AXPs/SGRs: Magnetars or Quark-stars? 
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Abstract 

The magnetar model and a solid quark star model for anomalous X-ray pulsars/soft 

■ gamma-ray repeaters (AXPs/SGRs) are discussed. Different manifestations of pulsar- 
like stars are speculated to be due to both their nature (e.g., mass and strain) and 
their nurture (ambience, and consequently the type of accretion) in the solid quark 

■ star scenario. Relevant arguments made by the author's group, including a debate 

X) 

o 



on solid cold quark matter, are briefly summarized too. 



1 Historical notes on magnetars and quark-stars 
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Since the discovery of radio pulsars about 40 years ago, more kinds of pulsar- 
like stars have been discovered, which are distinguished by their various man- 
ifestations, i.e., accretion-driven X-ray pulsars, X-ray bursters, AXPs/SGRs 
(anomalous X-ray pulsars/soft gamma-ray repeaters), CCOs (central com- 
pact objects in supernova remnants), DTNs (dim thermal neutron stars), 
and RRATS (rotating radio transient sources). Among these, the enigmatic 
sources, AXPs and SGRs, share three common surprising features: (i) all the 
sources have long spin periods (5 — 12 s) and persistent X-ray emission with 
luminosities well in excess of their spin-down powers, (ii) no clear evidence for 
them in binaries is obtained, and (iii) they show X-ray bursts, even giant flares 
saturating space detectors (the typical luminosity is greater than 10 6 times of 
the Eddington one). Four SGRs (4 confirmed, 3 candidates) and eight AXPs 
(8 confirmed, 1 candidate) have now been discovered (see the catalog updated 
online at '|http: / / www.physics.mcgill.ca/ ^pulsar/magnetar/main.html ' ) . 



Rapid rotation was generally thought to be the only e nergy source for pulsar 



emiss ion soon after the discovery of radio pulsars (e.g.. Manchester fc Tavlorl. 



19771 ) until the discovery of accretion-powered pulsars in X-ray binaries (jPringle fc Reed . 
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19721 ). However, AXPs/SGRs have long spin periods (thus low spindown power, 
their X-ray luminosities are much larger than their spindown powers) and no 
binary companions, which rules out spin and accretion in binary system as the 
power sources for the emission. It was then proposed that SGR-like bursts as 
well as the persistent X-ray emission could plausibly be the result of field de- 
cay of ultra-magnetic neutron stars if MHD dynamo action in the proto-stars 
i s very effective in cas e that the objects spin initia l ly at period s of ~ 1 ms 
( IDuncan fc Thompson! . Il992l ; iThompson fc Duncan! . Il993l . Il995l ) . Because of 
the starquakes in the crusts of normal neutron stars, a self-induction electric 
field is created. The strong electric field could initiate avalanches of pair cre- 
ation in the magnetosphe re and certainly accelerate particl es, resulting in a so- 
called magnetar corona ( Beloborodov fc Thompson! . 120071 ) . from which high- 



energy bursts could be observed. The power source of A XPs/SGRs is actually 
the m agnetic energy through field-reconnection there (IWoods fc Thompson! . 
20051 ) . This magnetar model is very popular nowadays in the astrophysical 
community. 



Is there any other kind of possible power source for AXPs/SGRs? This fun- 
damental question, however, could be related to a much more elementary and 
very important puzzle: What is the real nature of pulsar- like stars? (That is, 
are they really neutron st ars or quark stars? For a general review, see, e.g., 
Horvathl . 20061 ; Xu . 2006a ). In fact, alternative po wer sources wer e also pro- 



posed for AXPs/SGRs in a solid quark star model (jHorvathl . 120051 ; IXu et al. 
20061 ). Surely the properties of matter at high density but low temperature 
have been of great interests in both theoretical and observational aspects, es- 
pecially in the case of compact stars (see a discussion in §4). Quantum chromo- 
dynamics (QCD) is believed to be the underlying theory for the elementary 
strong interaction, but its strong coupling at low energy scales has troubled 
physicsists for a very long time. Lattice QCD simulations and effective QCD 
models are ways for us to study the non-perturbative effects, whereas lattice 
QCD seems not accessible for calculating supranuclear matter at low tempera- 
ture because of the still unsolved sign problem. Therefore, physicists can only 
approximate the state of matter at high density and low temperature with 
effective models. 



As for the inner constitution of pulsar-like stars, nuclear matter (related to 
neutron stars) is one of the speculations even from the Landau's time, but 
quark matter (related to quark stars) is an alternati v e due to the fact of Q CD's 
asymp totic freedom dBodmerl. Il97ll; IWittenl . Il984j ; iFarhi fc Jaffd . Il984f ) . Re- 
cently, iJaikumar. Reddy fc Steiner (2006) even proposed a crust comprised of 
nuggets of strange quark matter embedded in an uniform electron background, 
and found that a much reduced density gradient and negligible electric field 
would exist in such a quark star. It is a pity that we are now not able to de- 
termine confidently which state exists in reality due to the difficulty of QCD 
calculations. Therefore, we have to focus on astrophysical observations in the 



2 



research. Based on the Planck-like spectrum without atomic features and the 
precession properties of pulsar-like stars, a solid quark matter conjecture was 
suggested (jXu . 20031 ). Additionally, other features naturally explained within 
this model could possibly include sub-pulse drifting in radio emission, glitches, 
strong magnetic fields, birth after a successful core-collapse supernova, detec- 
tion of small bolometric radii, and transient bursts of GCRT J 1745— 3009 



flXu et all 1199a IZhou et all 12004! ; IXul . l2005al : lYue et all l2006al : IZhu fc Xu 



20061 ). Besides, the observational features of AXPs/SGRs may also reflect the 



nature of solid quark stars, that will be focus on in this paper. Power sources 
in solid quark stars include the elastic energy, the quake-released gr avitational 



energ y, and the accretion-released gravitational energy of fossil disk ( jXu et al. 



20061 ) . The hard X-ray bursts of AXPs/SGRs are proposed to be the results 
of AlSq (accretion-induced starquakes) in the model. 



2 Possible issues challenging the magnetar idea 



Conventional models for AXPs/SGRs are magnetars, a kind of neutron stars 
with surface magnetic fields in the range ~ (10 14 — 10 15 ) G, especially for 
understanding three supergiant flares (with released energy of ~ io 44-47 erg) 
observed in SGRs. The quiescent X-ray emission with luminosity ~ 10 34-36 
erg/s as well as the bursts of AXPs/SGRs are supposed to be powered by 
magnetic field decay. However, there are still some observational or theoret- 
ical arguments which may challenge the idea, although magnetars are really 
popular in the astrophysical community nowadays. 



In fact, it is still debatable to estimate the magnetic fields (oc V PP) of 
AXPs/SGRs simply from their period (P) and period derivative (P), since 
AXPs/SGRs could not be braked dominantly by the magnetodipole torque. 
AXPs/SGRs were suggested to be normal-field pulsar-like stars which are in 
an ac cretion propeller phase (lAlparl . l200ll ; IChatterjee. Hernquist &: Narayanl . 
20001 ) . But the difficulty in the latter view point is to reproduce the irregular 
bursts, even the superflares. That difficulty could be overcome if bur s ts are 
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proposed to be the results of starquakes of solid quark stars (lOwen 
Xu et all 120061 ) . In addition, some problems with the magnetar scenario may 
arise from recent relevant research, which are summarized as following. 



(i). The superstrong fields of magnetars are supposed to be created by MHD- 
dynamo action of rapid rotating protoneutron stars with spin period < 3 ms. 
The Poynting flux and the relativistic particle ejection of such a star should 
power effectively the supernova remnants. Numerical simulations of supernova 
remnant expansion with internal magnetars did shown a faster expansion when 
the energy injecte d into the supernova re mnant by magnetar spin-down is 
taken into account (lAllen fc Horvathl . 120041 ). Such energetic remnants are then 
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expec ted in magnetar models, but have not been detected (IVink fc Kuiperl . 
2006h . 



(ii). Dust emission around pulsar-like stars (e.g., AXPs) was proposed to test 
observationally th e propeller scenarios of quark stars with Spitzer or SCUBA 
jXil l2005aL l2006ch . Actually a recent discovery of mid-infrared emission from 
a cool disk around an isolated young anomalous X-ray pulsar, 4U 0142+61, 
has been reported (I Wang. Chakrabarty fc Kaplan! . 120061 ) although it is still a 
matter of debate whether significant propeller torque of fallback matter acts 
on the star. 



(iii). In order to understand the thermal Planck-like spectra in the conven- 
tional model of neutron star atmospheres for DTNs or CCOs, strong mag- 
netic fields (and possibly r apid rotat ion) are suggested. However, a problem 



inherent in these efforts is (jXul . 120041 ) : How to calm down the magneto spheric 



activities so that no AXP/SGR-like persistent and bursting emission could be 
radiated? This problem exists also for young "hi gh-field" radio pulsars since 
no magnetar-like emission is detected from them (IGonzalez et all 120061 ) . The 
facts could be understandable if both AXPs/SGRs and "high"-field pulsars are 
braked effectively by propeller torques of fossil disks, with the relative accre- 
tion rates (being indicated by the ratio of co-rotation radius to magnetospheric 
one) to be higher in the former than in the latter. 



iv). It is still a matter of debate whether the absorption feature s in SGR 1806- 
20 ca n be interpreted by proton or electron cyclotron resonance (jXu. Wang fc Qiao . 
20031 ). The field is only ~5x 10 11 G in the context of electron cyclotron lines. 



(v). The pressure should be very anisotropic in a relativistic degenerate neu- 
tron gas in equilibrium with a background of electrons and protons when the 
magnetic field is stronger than the critical field (~ 10 13 G). The vanishing of 
the equatorial pressure of the gas might r esult in a transverse collapse, an d 
a stable magnetar could then be unlikely ([Martinez. Rojas fc Cuestal . 120031 ). 
Such kind of quantum-magnetically induced collapse could also be possible 



in qu ark stars with ultra strong mag netic fields (> 10~ 19 G, Martinez et all 



20051 ). In this sense, quark matter would tolerate much higher magnetic field 
than neutron matter, and no neutron star with 10~( 14-15 ) (i.e., no magnetar) 
could form stably. Nevertheless, more theoretical investigations of the equilib- 
rium of nuclear or quark matter in general relativity are necessary. 



In addition, the implied transverse velocity of an AXP, XTE J1810-197, in- 
dicates that its kick velocity is completely o rdina ry (< 400 km/s) though 
the AXP manifests strikingly (IHelfand et all 120071 ). It might be not easy to 
reproduce normal kicks for neutron stars with very fast rotation (initial pe- 
riod P < 3 ms) and much high magnetic fields (B ^> 10 13 G) in the magnetar 
model, since the initial magnetic field and the initial spin are generally consid- 
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ered to be key parameters for kicks (jLail . 12004 ). Considering these criticisms, 
we think that alternative ideas for understanding AXPs/SGRs are welcome. 
We note that the quark-star model could be competitive, in which both grav- 
itational and elastic energies would power the photon and particle radiation 
of AXPs/SGRs. 



3 AXPs/SGRs: quark-stars with fossil disks? 



Why do pulsar-like stars behave so diversely? As noted by Kul (j2006d ). two 
factors could be very important for their different manifestations: environ- 
ment and stellar mass. A dense environment may be responsible for leaving 
a remnant-nebula after a supernova and for resulting possibly in a fallback 
accretion (e.g., via disk) around the compact star, whereas starquakes could 
more probably occur in massive compact stars. Bursts, as the results of quakes, 
could frequently appear in an accreting solid quark star (even low-mass) be- 
cause of AlSq (accretion-induced starquakes). In this context, various mani- 
festations of pulsar-like stars would be speculated to be related to quark stars 
with different nature and nurture (fossil disks , quakes, and m asses; see Table 
1). According to the calculation of Fig. 2 in Ku et al.l (120061 ). more gravita- 
tional energy could be released during quakes of solid quark stars with higher 
masses, and we may expect that SGRs emitting giant flares could be quark 
stars with ~ (1 — 2)M Q . 



By assuming that AXPs/ SGRs and "high"-fi eld radio pulsars could be torqued 
effectively by fossil disks, lYue et all fl2006bh suggest that SGR 1900+14 and 
AXP IE 1048.1-5937 might be quark star with ~ 1O" (1 - 2) M mass if all AXPs 
and SGRs are likely grouped together in their Fig. 1. One of the "high"-field 
pulsar, PSR J1846-0258 in the supernova remnant Kes 75, has only been 
detected in the X-ray band, with a much higher ratio of X-ray luminosity 
to the spindown energy loss than that found in other young radio pulsars 
(e.g., the Crab pulsar). The timing of PS R J1846-0258 behaves lik e that of 
radio pulsars, with a braking index of 2.65 ([Livingstone et al.l . 120061 ). Though 
the emission luminosity is smaller than the spindown power, it is still possible 
that PSR J1846-0258 would be actually accretion-powered if the real magnetic 
moment is very low. In this case, the monopole-induced potentials drop in the 
open field region of PSR J1846-0258 might not be high enough (due to a 
much lower real magnetic field) to produce significant pair plasma for radio 
emission. Also the potential drop of AXPs/SGRs could be too low to initiate 
avalanches of e ± pairs in their quiescent states. However, it is worth noting 
that radio activity would still be possible for AXPs /SGRs after outburs t, as in 



2006|) 



since 



the case of radio emission from XTE J1810— 197 (jCamilo et al.l . 
strong pair plasma could form due to quake-induced E\\ (electric field parallel 
to magnetic field). In a conclusion, AXPs and "high" -field radio pulsars might 
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Table 1 

Speculated characters of pulsar-like stars in the quark star scenario. 



manifestation 



propeller torque 



starquake 



mass 



Radio 
pulsars 


normal 
millisecond 
part-time 
"high" -field 


negligible 
negligible 

THCGl i~\ 1 G u ( 

lUooll CLloK I 

significant 


glitches 
microglitch* 

glitches 


some low? 
most low? 

1UW ICljlll^Mi 


Accretion- 
powered 


X-ray pulsars 
X-ray bursts 


strong accretion 
strong accretion 


a frequency glitch' 




CCOs 


significant 


no-detection 


low? 


DTNs 


weak accretion 


no-detection 


low? 


AXPs 


substantial 


AlSq*, frequent 1 ' 




SGRs 


substantial 


AlSq, more frequent 


~ Mq for those 
with superflares 



* Only one microglitch with 10 11 c hange in fractional f reque ncy was detected in a 
millisecond pulsar, PSR B1821-24 (ICognard & Backer! , liool ) . 

t Still one glitch in an accreting pulsar, KS 1947+300, has been reported 
during outbursts, with a fractional change in frequency of 3.7 x 10~ 5 to be 
significantly larger than the valu es observed in the glitches of radio pulsars 
(jGalloway. Morgan & Levinel . l2004l ). 

* Accretion- induced starquakes (AlSq). 

^ Large glitches are detected from two AXPs, nearly simultaneously with burst: IE 
2259+586 ( Aylv ~ 4 x 10~ 6 ) and CXOU J164710.2-455216 (Au/u = 6.5 x 10~ 5 ) 
(jlsrael et all I2007I ). 



be in a same group, with different magnetic moment but similar appearances. 



The magnetospheric electrodynamics of isolated rotation-powered pulsars is 
studied extensively in the literatures, including pair production, particle ac- 
celeration, and E 1 1 -solutions. However, the corresponding dynamics of a rotat- 
ing quark star in a propeller phase, with significant accretion onto the polar 
cap, could be in a very differen t scen ario. Ruderman-Sutherland-type vacuum 
gap (jRuderman &: Sutherland! Il975l ) may not exist (thus no sparking occurs) 
above polar caps since the hot polar cap bombarded by accreted nucleons 
could emit sufficient charged mesons and/or hadrons in the QCD-phase con- 
versio n process. A space charge-limited flow model (lArons fc Scharlemannl . 
19791 ) may adaptable in this case, where the mesons and hadrons are accel- 
erated effectively by En and may annihilate finally (or do by other interac- 
tions) into ultra-high energy photons and neutrinos in magnet ospheres. TeV- 
photon emission observed could be understood in this picture, which h a s also 



been noted by many authors (e.g., IZhang et all 120031 ; iLink fc Burgiol . 12005 



Neronov fc Chernyakoval . l2006l ). Certainly, the related dynamics needs further 
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investigations in order to study the detail astrophysics there. 



We note that fossil disk could affect the spindown torque in two aspects, (i) 
Propeller torque attributes to the MHD-coupling between accretion flow and 
co-rotating particles at the magnetospheric radius (r m ). (ii) Particles (acceler- 
ated by E\\) with higher number density would be ejected (and the star spins 
down faster) if r m is smaller than the light cylinder radius (r\ c ; the open field 
line region is then determined by r m rather by n c ). In this case, the potential 
drop in the open field lines is higher than that expected for pure rotation- 
powered pulsars. 



Conclusive evidence for quark-stars ? \X.xa ( I2006d j had addressed that following 
observations could test for the quark star (even with low mass) idea: (i) to 
detect dust emission around pulsar-like stars; (ii) to geometrically determine 
the radii of distant pulsar-like stars by advanced X-ray or UV-band facilities; 
(iii) to detect gravitational waves from pulsar-like stars; and (iv) to search sub- 
millisecond pulsars (especially with periods < 0.5 ms). Although it should cer- 
tainly be very clear evidence for quark stars if we could detect sub- mill i secon d 
pulsars by advanced radio telescope (e.g., the FAST) in the future (jXul . 120071 ) . 
we note that to detect dust emission form pulsar-like stars is another way to 
find such evidence. The maximum elastic quadrupole deformations sustain- 
able by normal neutron stars and sol id strange q uark stars might result in 
distinguishable maximum ellipticities (lOwenl . 120051 ). and the gravitational ra- 
diation may show the nature of pulsar-like stars. Though we shouldn't detect 
gravitational wave because of th e low spindow n powers in the second LIGO 
science run ( lAbbott et all 120051 ; lOwenl . 120061 ). a result of search for known 



pulsars in future LIGO run might only show upper limits of masses (or radii) 
for the targe ts since the gravitational wave radiation could be mass-dependent 
jXul l2006bh . It is also worth noting that gravitational wave emission associ- 
ated with the e nerge tic superflares of S GRs may be gathered from the LIGO 
data ( Horvathl . 2005 ). and that, in fact, LIGO Scientific Collaboration ( 2007 ) 
had put an upper limit on the gravitational wave emission, 7.7 x 10 46 erg, for 
the 92.5 Hz QPO of SGR 1806-20. In addition, it is still not sure whether the 
small radii derived from Planckian thermal radiation indicate the real stellar 
surfaces or only the hot polar caps, and it is then very necessary and important 
to observe confidently the emission from global surface of pulsar-like stars in, 
e.g., soft X-ray or UV bands. To find a binary system of low mass quark star 
(formed via AIC) and white dwarf could also be expected in the quark star 
model, which might be an effective way to search quark stars with very low 
masses. 
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4 Conclusions and Discussions 



One of the daunting challenges nowadays is to understand the fundamental 
strong interaction between quarks, especially in the low-energy limit, since the 
coupling is asymptotically free in the limit of high-energy. Though quark mat- 
ter is predicted in QCD because of asymptotic freedom, astrophysical quark 
stars would be still in a low-energy region where quarks are coupled non- 
linearly. In this sense, AXPs/SGRs (including other pulsar-like stars) would 
be ideal laboratories for us to study the non-perturbative effects of QCD. Phe- 
nomenologically, a solid state of cold quark matter was suggested, on which 
understanding the peculiar features of AXPs/SGRs could be based. It is sug- 
gested in this paper that the various manifestations of pulsar-like stars could 
imply different environments and inner structure of the sources, and that the 
bursts of AXPs/SGRs could be the results of accretion-induced starquakes 
(AlSq) of quark stars in their propeller phases. A solid quark star quakes as 
stress energy develops due to spin-down, or cooling, or even general relativistic 
effect. 



There exist other versions (e.g.. lOuyed. Elgaroy &: Dahlel . l2006al ; IOuyed. Leahy fc Niebergall . 
2006al lbh of quark star models for AXPs/SGRs, which may show the neces- 
sary of introducing astrophysical quark ma tter for understanding this kind o f 
sources. Assuming a superconducting state. iNiebergal. Ouyed fc Leahyl (120061 ) 
calculated the magnetic field decay due to the expulsion of spin-induced vor- 
tices, and found that the magnetic field strengths and periods remain almost 
unchanged for typical parameters. Though glitches in the model of Ouyed et 
al. could be explained with intricate details, glitches and quasi-periodic oscil- 
lations (QPOs) during superflares mig ht hint that the cold quark matter could 



be in a solid state. A simple toy-model (iGlampedakis. Samuelsson fc Andersson 
20061 ). with the inclusion of global MHD-coupling between the elastic crust 



and the fluid core of a normal neutron star, was suggested to explain QPOs of 
SGRs, but a similar model in the solid quark star scenario is worth to make 
in order to understand the detail observations and to know observationally 
the features of solid quark matter (e.g., the shear modulus), since MHD waves 
in the global magnetospheres could be excited by oscillations of quark stars 
after their quakes. It is also very necessary to calculate the quake-induced E\\ 
(electric field parallel to magnetic field) and particle acceleration in both open 
and closed field line regions in order to model the magnetospheric activity of 
AXPs/SGRs after a s tarquake. Two large glitches (Av/v ~ (10" 5 - 10~ 6 )) 



are detected in AXPs (jlsrael et all 120071 ; iMuno et all 120071 ). and it would be 



very urgent to search glitches in SGRs since larger glitches could be natural 
in this AlSq scenario. 



Actually, from a viewpoint of condensed matter physics, symmetries break as 
energy-scale decreases so that different states of matter are in the Universe. 
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Cold quark matter was supposed to be approximated by ideal Fermi gas when 
the interaction between quarks is neglected for matter to be asymptotically 
dense. But a superfTuid state, in which the global U(l) symmetry is broken, was 
also proposed since any attractive interaction between quarks may render the 
Fermi surface unstable. It is further suggested that a crystalline structure (i.e., 
the translation symmetry is broken) may form in this color superconductivity 
quark matter, and that the shear modulus could be 20 to 1000 time s larger 



than those of neutron star crusts (jMannarelli. Rajagopal &: Sharmal . 120071 ) 



However, another possibility of condensation in position space could still not 
be rule out. Quark clusters would form in this scenario, and quark matter 
would be solidified at low temperature. An investigation on the initial spin- 
evolution may show signatures of quark star solidification. Star's cooling and 
gravitational radiation would affect its viscosity and rotation. 

From sQGP to solid quark matter? The non-linear complexity of strong in- 
teraction at low energy could lead to two hot points in recent researches: the 
physics of relativistic heavy ion collision and the study of pulsar-like star's 
structure. A key point in the quark-star model suggested for AXPs/SGRs is 
the conjecture that quarks would be clustered in strongly coupled quark-gluon 
plasma (sQGP), especia lly in the supranuclear matter with low temperature 
( Xu , 2003 . 2004 . 2005al Jbl) when the interaction parameter T (defined by the 



ratio of the potential energy to the kinetic energy) is very high. Solid quark 
matter could be very likely in case that the thermal kinematic energy of quark 
clusters is much lower than the interaction energy between the clusters (i.e., 
r ^> 1). In fact, it is recently realized T- value of the QGP at RHIC should 
be much highe r than that predicted by the perturbative theory of QCD (e.g., 
Shurvad . l2006h . Therefore, the QGP created at RHIC is also sQGP, and one 



may expect an even higher T for cold quark matter (i.e., temperature to be 
much lower than 1 MeV). Glass-like or crystal-like solid quark matter would 
then be possible in case of keV-temperature. One way to test the phenomeno- 
logical interaction-models for sQGP at laboratory could be calculating the 
quark clusters by qMD (quark molecular dynamics), where the colored and 
flavored quarks c ould be treated as sem i-classical particles interacting via a 



Cornell-potential (IHofmann et al.l . 120061 ) 



V ( r ) = h nr. 

4 r 

Strong interaction via this potential may favor formation of quark clusters, al- 
though the interacting potential between quark clusters could be in a different 
type. The experimental yield of multi-quark clusters (pentaquark, teraquark) 
may improve our underst anding of sQGP. Rece nt achievements about classi- 
cal sQGP are presented in I German et al. fl2006al ]bh. Therefore, we expect that 



the idea of quark clustering in cold quark matter would probably be tested 
experimentally. 
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The beauty of non-perturbative QCD? The uncertainty of determining the 
composition of pulsar-like stars depends actually and strongly on the non- 
perturbative nature of QCD at low-energy level, which is unfortunately not 
well understood. This is the real reason that we cannot make sure pulsars 
are neutron or quark stars according only to today's QCD calculations. Pes- 
simistically, it might be impossible forever to know the state of supranuclear 
matter by pure QCD calculations; and we should then have to recognize the 
importance of observations and phenomenology in studying pulsar's inner 
structures. Investigations in such a way may discover the natural "beauty" 
of non-perturbative QCD. As in another but simple kind of non-perturbative 
dynamics, the non-linear phenomena in fluid dynamics (such as turbulence 
and chaos), which are surely not able to be determined by the famous Navier- 
Stokes equation, do exist in nature and are beautiful. In this sense, it is not 
against the rule of nature's beauty to study observationally and phenomeno- 
logically the elementary strong interaction at low energy. 
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